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We present a focused study of a predictive unified model whose measurable consequences are 
immediately relevant to early discovery prospects of supersymmetry at the LHC. ATLAS and CMS 



! have released their analysis with 35 pb ^ of data and the model class we discuss is consistent with 

Ph. 

^1 this data. It is shown that with an increase in luminosity the LSP dark matter mass and the gluino 

D ■ 

r~| ' mass can be inferred from simple observables such as kinematic edges in leptonic channels and peak 



values in effective mass distributions. Specifically, we consider cases in which the neutralino is of 
low mass and where the relic density consistent with WMAP observations arises via the exchange 



> 
OO 

' of Higgs bosons in unified supergravity models. The magnitudes of the gaugino masses are sharply 



limited to focused regions of the parameter space, and in particular the dark matter mass lies in 
the range ~ (50 — 65) GeV with an upper bound on the gluino mass of 575 GeV, with a typical 
mass of 450 GeV. We find that all model points in this paradigm are discoverable at the LHC at 
^ ■ y/s = 7 TeV. We determine lower bounds on the entire sparticle spectrum in this model based 

. on existing experimental constraints. In addition, we find the spin-independent cross section for 

neutralino scattering on nucleons to be generally in the range of a% = 10-46±i 

cm with much 

higher cross sections also possible. Thus direct detection experiments such as CDMS and XENON 
already constrain some of the allowed parameter space of the low mass gaugino models and further 
data will provide important cross-checks of the model assumptions in the near future. 
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I. INTRODUCTION 

Unified models of supergravity witli gravity mediated breaking of supersymmetry |l| 
extend tlie standard model of particle physics and are being tested with the Large Hadron 
Collider experiments at CERN. As a consequence of the breaking of supersymmetry, one 
obtains soft masses and couplings of the form [l, 2| 

mi/2 = M3(Au) = M2(Au) = Mi(Au), (1) 



m, 



= <(Au) = m|(Au) = m^^^(Au), (2) 



Ao = A,„t,bA^v) , (3) 

where at the unification scale, Ay ~ 2 x 10^^ GeV, there are universal mass terms for 
the gauginos of SU{3), SU{2),U{1), denoted by mi/2, and universal mass squared terms 
for scalar fields denoted by rriQ (where Q (L) stands for squarks (sleptons)), and universal 
cubic (trilinear) couplings Aq which multiply the Yukawa couplings of matter fields to the 
Higgs fields. In addition, a (bilinear) soft Higgs mixing term proportional to fiQ of the form 
Bqi1q[HiH2 + h.c.) arises from the superpotential, where H2{Hi) are the Higgs doublets 
which give mass to the up quarks (down quarks and charged leptons). The constraints 
of electroweak symmetry breaking allow the determination of \fi\ (where fi is /io at the 
electroweak scale) in terms of Mz and further one makes the replacement of Bq by the ratio 
of the Higgs vacuum expectation values tan /3 = {HS)/{H^) leaving minimally 4 parameters 



and one sign needed as input to define the model 



mo, mi/2, Aq, tan/3, sign(/i) . (4) 

Through renormalization group evolution, one computes the predictions for all the masses 
of the superpartners and their couplings to each other and to the standard model fields^. 

Models of supergravity address fundamental questions in particle physics, such as the 
gauge hierarchy problem, the breaking of electroweak symmetry, and the unification of 
strong and electroweak forces. In addition, such models also provide a compelling dark 
matter candidate; the lightest supersymmetric particle (LSP). In particular, the neutralino 
is a linear combination of gauginos and Higgsinos as follows: 

= niiB + ni2W + riisHi + nuH2, (5) 



For recent reviews see: 



where B is the bino, W is the wino and H12 are the Higgsinos. The neutrahno can have 
the right cross section and mass to provide a natural candidate for the observed density of 



cold dark matter (CDM) in the universe. According to the analysis in j?], the latter has the 
value 

f]cDM/i^ = 0.1120 ±0.0056 . (6) 

Here h is the Hubble constant, Hq, in units of 100 km/s/Mpc, and under the assumption 
that ficDM = ^x?' '^^^ ^x? ~ Pxll where the neutralino density p^o is in units of the 
critical density pc = SHq/ISttG) ~ 2 x 10^^^ h"^ g/cm^. The measurement of the relic density 
together with a variety of results from collider experiments provide strong constraints on 
models of new physics. 

In this paper we study a particular region of the unified supersymmetric parameter space 
which satisfies all the existing experimental and astrophysical bounds and is testable in the 
very near future. We focus on the region where the neutralino has a mass in the range 
~ (50 — 65) GeV. In this mass range, which is above the Z-pole, when 2m^o < m/j, in those 
models that are unconstrained by present experimental data, the relic density of neutralinos 
is largely governed by the presence of the light CP even Higgs pole (h-pole) js, 9| through 
annihilations in the early universe, schematically: 

XiXi ^ h^bb, Tf, cc... i2m^o < nih) (7) 

arising from the resonance; however, other channels can contribute in general. Additionally, 
when 2m^o > rrih the relic density can also be achieved via 

^h,H,A-^ ff (8) 



through the s-channel where the heavier Higgses can play the dominant role [8|. Such 
annihilations can lead to effects on the relic density when the mass of the pseudoscalar 
is light, of order a few hundred GeV, which corresponds to the case of large tan/3. Our 
analysis will find results consistent with a large range of tan /3 ~ (3, 60) with the possibility 
of both a heavy and a light pseudoscalar. We will refer to the collective region of the 
parameter space, with |m^o — m/i/2|max ^ 0{5) GeV as the "Higgs-pole region". 

With universal boundary conditions at the unification scale, the mass range of the neu- 
tralino is confined by mass limits on the other particles in the spectrum. In particular the 



light chargino has a bound from LEP of m^± > 103.5 GeV [l^. It is known that in models 



with the minimal supersymmetric field content the light CP-even Higgs mass has an upper 

n 

bound of roughly nih ^130 GeV [11]. The Higgs mass is bounded from below by direct 



searches at LEP [l^ and, more recently, at the Tevatron [13|. We will use a conservative 
lower bound of ruh > 110 GeV to allow for the theoretical uncertainty in computing the loop 
corrections to the Higgs mass. We note that a stricter imposition of nih > 114 GeV would 
narrow the space of models but has little impact on our generic conclusions. Specifically, 
the low mass gaugino models we study in the Higgs-pole region will correspond to light 
neutralino dark matter in the range 

52 GeV < m^o < 67 GeV (9) 

that yields the correct relic density and obeys all other experimental constraints subject to 
the boundary conditions of Eq. . 

Here we will show explicitly with a dedicated study that this class of low mass gaugino 
models should either be found or ruled out with early LHC data if the expected luminosity 
of ~ few fb~^ is reached at ^/s = 7 TeV. In addition, we will discuss current and upcoming 
dark matter direct detection experiments which also have the possibility of detecting the 
neutralino LSP in these models. 

The reason the models in the Higgs-pole region can be tested soon is that several impor- 
tant mass scales are low enough to be within the discoverable reach of LHC-7. It is known 
that in minimal supergravity models the following scaling relation amongst the neutralino 
LSP, the chargino, next to lightest neutralino, and the gluino masses are satisfied ^ 

2m^o ~ m-± ~ m^o ~ -rria ■ (10) 

Al Xi A2 ^ a ' 



For a precise determination o: 



to the gaugino masses 



the scaling relations above one must include loop corrections 
2l[ |. Eq. flTU]) typically holds for a very pure bino LSP; whereas 
the scaling relations receive significant corrections when the LSP eigenstate has a non- 
negligible Higgsino component. The constraint of Eq. ( ITO|) . which we will generalize, is an 
important guide regarding the types of signatures at the LHC for this class of models. In 
what follows we will take the scaling assumption to mean that the mass relations of Eq. ( !T0|) 
(or the generalization thereof, which is included in Eq. ( [T3l) in what follows) hold to a good 
approximation. 



This relation holds for the case when /k^ ^ Mf, Mf, M| all taken at the electroweak scale. 
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Remarkably, in the literature there are rather few studies of the impact on LHC physics 
from this Higgs-pole region with correspondingly low mass gauginos; only recently has it 
seen some attention. Thus, some aspects of the minimal supergravity models where the relics 



annihilate near the light CP-even Higgs pole have been discussed in Ref. 



14- 



under the mass hierarchy denoted by mSP4 (supergravity mass pattern 4) IJ, |l5|, where. 



li^hich fall 



in particular, a clean edge in the dilepton invariant mass in this model class was noted 



in Ref. 15|]. In addition, the very recent work of Ref. 19| studies electroweak symmetry 
breaking in an overlapping class of models with a focus on the /i parameter and radiative 
breaking. 



Some of our observations and emphasis here have overlap with Refs. [171] and some 



are rather different. In Ref. 17| emphasis was given to explaining the CDMS II results 
and predictions for the XENON data, and in doing so, a slice of the parameter space was 
studied where tan /3 = 50 and Aq was fixed for a few choice values, while the analysis allowed 
for flavor violation, and thus constraints from 6 — )■ 57 and Bg — )■ fi~^^~ were not imposed. 
Our present analysis imposes these constraints and opens up new parameter space where 
all direct and indirect constraints are satisfied, and where the spin independent scattering 
cross section can lead to event rates that can be observed in the XENON detector. 

When all direct search limits and indirect constraints on the parameter space are imposed 
a number of robust mass relations are predicted. The main points emphasized in this work 
are as follows: 

1. Two key observables which are directly measureable at the LHC: the peak in the 
effective mass distribution as well as the dilepton invariant mass edge are shown to 
be strongly correlated in these models. A first determination of the gluino mass can 
he measured from the peak value of the effective mass distribution and the dark matter 
mass can simultaneously he inferred from the dilepton edge due to the predicted scaling 
relations in the gaugino sector given in Eq. (Q] 



2. The recent CMS and ATLAS data with 35 pb^^ of integrated luminosity [26|, |27| do 
not yet provide constraints on the models discussed in this paper. In the Higgs-pole 
region, even though the gluino has a low mass, the 2'^'^ generation squark masses are 
larger than 1 TeV and typically of order several TeV which is the main reason these 
models remain unconstrained by the CMS and ATLAS data (the gluino mass bounds 
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in the recent ATLAS analysis do not apply to our models). However, we will show 
that with increased luminosity they will begin to probe such models. 

3. The gluino has a low mass which is tightly constrained to lie in the range 400 GeV ^ 
rrig ^ 575 GeV, with most points having^ rrig ~ 450 ± 20 GeV. The mass splitting 
between the gluino and the lighter gauginos is appreciable. Thus should this model 
class be realized in nature, the production of jets from the gluino must be seen at the 



LHC yi= 7TeV with about 1 fb^^ of data [16j,|22|,|23j,y, l25j,(l9|. 



4. The chargino mass is bounded from below by the LEP search limits and from above 
by theory, m^± ^ 130 GeV, with the second heaviest neutralino being effectively 
degenerate with the lightest chargino. This suggests that the associated production 
of X1X2 is sizeable and may reveal itself in multilepton channels, in particular the 



trilepton (3L) channel 28|, |29|]. The large SUSY breaking scalar masses in the models 



imply that the current bounds from the Tevatron do not yet constrain the models. 

5. There is a sizable region of the parameter space in which tan /3 can be large and 
the pseudoscalar Higgs boson is relatively light. Such model points may allow for 
simultaneous reconstruction of nig and in early LHC data collection. 



6. The constraints from the CDMS and XENON data [30|, l3l|] on the spin independent 
scattering cross section of neutralinos on nucleons is complimentary to searches for 
the CP-odd Higgs at the Tevatron and at the LHC. In fact, for some models in the 
parameter space the XENON data already constrains models that will be tested in 2011 
and 2012 at the LHC. We find many candidate models that yield large event rates in 
upcoming dark matter direct detection experiments. 

As an aside, we note that the neutralino annihilation rate we consider is too low to 
produce observable cosmic signatures of positrons, antiprotons, or gamma rays; hence recent 
experimental bounds from a variety of cosmic ray experiments are not a concern. In principle 
one could boost the annihilation cross section in a number of ways in order to reach the 
sensitivity of the experiments, but that approach is not considered here. 



This is the Gaussian peak (i.e. mean) and Gaussian width (i.e. 1 standard deviation) 
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Thus in this work we study a dense region of the parameter space of minimal supergravity 
models where the LSPs have low mass that also have low mass gluinos which will be tested 
at the LHC in the very near future. In addition, we find a bound on the Higgs sector from 
the XENON data. We explore the connection between these models and what the LHC, 
the Tevatron, and the dark matter scattering experiments can observe. The prominent 
signatures of the models under full collider simulation are discussed in detail in what follows. 



II. ANALYSIS OF THE PARAMETER SPACE AND SPARTICLE MASSES 



In this section we describe our targeted parameter scan over the minimal supergravity 
parameter space for the low mass gaugino models that lie in the Higgs-pole region. We will 
illustrate the various constraints we have imposed on the models, from astrophysical relic 
density as well as accelerator bounds. From the results of our survey of parameter space, 
we then obtain the viable range for sparticle masses and the relations between them. 

In the analysis that follows we compute the thermal relic density as implemented in 
MicrOMEGAs 2.4 32|. We demand that the resulting value of the cold dark matter relic 



density flcBMh"^ = f^yo/i^ satisfy 



xr 

2 



0.08 < n^oh' < 0.14 . (11) 

The spread in f|TT]) around the WMAP band is chosen to allow theoretical uncertainties 
and sensitivity to the top pole mass, both of which enter in the sparticle spectrum under 
renormalization group flow and radiative electroweak symmetry breaking. 

Our targeted parameter scan over the minimal supergravity parameter space is described 
in what follows. For models in which the gaugino masses are given by a universal parameter 
mi/2 at the scale Aqut — 2 x 10^^ GeV, the analysis of Ref. {2!, Q found that Eq. ( 1T0|) is 
consistent with m-± ~ m^o ~ (0.9 ± 0.1)mi/2; thus in the interest of obtaining models with 
low mass gauginos, we restrict mi/2 to the range 100 GeV < mi/2 < 175 GeV. The universal 
scalar mass was allowed to vary in the range 0.1 TeV < mo < 10 TeV with the upper bound 
representing a naturalness requirement on the models. The entire allowed range of tan (3 
was explored and the universal trilinear parameter Aq was allowed to vary over the range 
—4 < Ao/mo < 4. Throughout we take /i > and m^°^ = 173.1 GeV. Renormalization 
group evolution and calculation of the physical masses of the sparticles was performed using 
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SuSpect [3J] and SUSY-HIT [35| was used in the computation of branching ratios of the 
superpartners. 

Our survey resulted in 12,000 parameter sets, each defining a single model. All model 
points were required to satisfy the requirements of radiative electroweak symmetry break- 
ing. Accelerator constraints were applied as well. The most important bounds include the 
imposition of the higgs mass bound discussed in the previous section, and the bound on 
the chargino mass from direct searches for sparticles m-± > 103.5 GeV from LEP llOl. In 
addition a number of indirect experimental constraints were imposed, which include those 
from the Tevatron, Belle/BaBar/Cleo and Brookhaven experiments. Specifically we im- 



pose the conservative constraints (—11.4 x 
Br {Bs < 4.7 X 10-^ (90 % C.L.) 



36 



0-1°) <S{g^-2)< (9.4 X lO^^), see 
and 2.77 < Br (6 ^ 57) x 10^ < 4.27 



37|, 
39|. 



The indirect constraints were calculated using Mi cr Omegas, with the Standard Mode^ 



con- 



tribution in the last observable corrected according to the work of Misiak et al. 37|, |40 |. 
Finally, we require that the relic density satisfy Eq. f lTTj) . 
The models surveyed are consistent with 



Al 



mh/2 



< 7GeV. 



(12) 



with most models satisfying |m^o — mij2\ < 4 GeV. Therefore, post facto, Eq. f lTT]) and 
Eq. (fT2l) together provide an effective definition of what constitutes the Higgs-pole region. 
From this ensemble of models we find the mass relations 



nT'h = C(hm^o, 1.78 < ah < 2.25 
m-± = a-±m^o, 1.65 < a-± < 2.07 



X.2 X.2 Al ' 



1.70 < a.o < 2.07 

— X2 — 



rrig = agm^o, 7.34 < Og < 9.25 



(13) 



and the qualitative scaling relations in Eq. ( JTOj) can be replaced by the more quantitative 
relations 



rrih 



0.92 < /3-± < 1.17, 0.22 </3g< 0.29 . 



(14) 



The distribution of gluino masses for the models is well approximated by a Gaussian with 
a remarkably small width. In Figure [1] we plot the distribution in the dimensionless ratio 
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nig /ni^o 



FIG. 1: Distribution of the ratio ag = rrig/m^o from Eq. ()13p . The distribution is well approximated 
by a Gaussian characterized by ag = 7.86 it 0.209. The corresponding spread of gluino masses for 
the models simulated was found to be rrig = (451 it 19.5) GeV (quoted are mean values and one 
standard deviation about the mean). 

Ug = nig/m^o from Eq. f fT3|) . We see that in general the models produce a gluino mass of 

m-g = (451 ± 19.5) GeV (1 a). (15) 

Thus consistent with Eq. (11 01) one finds 

m-g/m^o = 7.86 ± 0.209 (1 a). (16) 

In Table [H we expand on the general ranges given in Eq. (1131) . For example, whereas in the 
previous paragraph and in Figure [1] the la error bars are quoted for the gluino mass, the 
full range of all gluino masses obtained in our survey is 

396 GeV < rUg < 575 GeV . (17) 

The upper bound for the gluino mass, consistent with a low mass neutralino has very 
important consequences for LHC searches as discussed in the next section. Another result 
of our analysis is that while the LSP is dominantly bino-like it can also have a significant 
Higgsino component as seen from Table [H 
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Predictions for the Sparticle Masses and LSP Eigencontent 



Mass Predictions (GeV) 


Eigencontent of the LSP 


110 < ruh 


< 126 


0.888 < mi < 0.996 




52 < m^o 


< 67 


-0.163 < ni2 < -0.016 




104 < m-± 

Xi 


< 131 


0.019 < ni3 < 0.396 


{Hi) 


396 < nig 


< 575 


-0.167 < ni4 < -0.006 





TABLE I: General predictions for the sparticle masses for the models with mo < 10 TeV satisfying 
all phenomenological constraints discussed in the text. It is further found that mg > 1.05 TeV, 
and the scalar masses are bounded as : m^^ > 323 GeV, m^^ > 706 GeV, nif^ > 484 GeV, rriq > 
1070 GeV, > 1050 GeV, and niA > 187 GeV. 



For the small values of mi/2 that lead to a light gaugino sector it is necessary to require 
large mo and/or tan /3 to satisfy the direct search limits on the light CP-even Higgs mass 
h. We therefore found that tan /3 ranges from about 3 to 60 and that typically mo is much 
larger than mi/2- Indeed in our survey an empirical lower bound of mo > 1.05 TeV was 
obtained. A large fraction of the models thus lie on the hyperbolic branch/focus point 
region 4lJ in which scalars are in the TeV range and /i is typically small. Consequently all 
the first and the second generation squarks and sleptons are significantly heavier than the 
gluino. In particular one finds the lower bounds rriq > 1070 GeV and m^ > 1050 GeV on 
squarks and sleptons of the first two generations. Third generation squarks and sleptons 
are also found to be generally heavy, though lower masses occasionally arise for certain 
combinations of Ao/mo and tan/3. Specifically we find the following lower bounds on third 
generation scalars: m^-^ > 323 GeV, m^^ > 706 GeV and mf^ > 483 GeV. 



We further note that the parameter for most of the models lies in the range 
300 GeV ^ /U ^ 700 GeV, though larger values are possible. The models with low fi can 
lead to a CP-odd Higgs mass rriA that can be quite light - particularly when the value of 
tan/3 is simultaneously large. We find a lower limit of m^ > 187 GeV over the ensemble 
of models studied. As we will see below, inclusion of the limits on the neutralino-proton 
spin independent cross section, o"|op, from the CDMS and XENON experiments further 
constrain the models. We discuss this in some detail in Section [IVl 
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Finally, one might ask if charginos with masses in the range 104 GeV < m-± < 131 GeV 
are already ruled out by direct searches at the Tevatron, given the recently quoted lower 
bounds of m-± ^ 150 GeV derived from the absence of trilepton events with large missing 



transverse energy 42h44||. Such a lower bound is due to the assumption of light slepton 
masses. However, as discussed above, the low mass gaugino models in the Higgs-pole region 
single-out scenarios in which the sleptons are generally very heavy, as in the "large mo" 



models analyzed by D0 



431]. Using Prospino2 45|] to calculate the next-to-leading order 



(NLO) production cross sections for the Tevatron at = 1.96 TeV we find, before cuts and 
efficiency factors, 

1.33 X 10-2 pb < a{pp ^ X2X^ )NLoBr(x^ ^ /^z/x?)Br(x^ ^ l+rxl) < 5.98 x 10"^ pb (18) 

after simply summing over all three generations of leptonic decay products, which is the 
maximal case, and this result is below the reported limits from the Tevatron 421-144 1 . 



III. SIGNATURES OF THE LOW MASS GAUGINO MODELS IN THE HIGGS- 
POLE REGION AT THE LHC 



To study the signatures of the low mass gaugino models at LHC-7 we simulate events 
at a/s = 7 TeV for a sample of 700 model points from the larger set discussed in the 
previous section. The standard model (SM) backgrounds considered were those used 



m 



24 



25| which compare well to those given in 23|. The SM b ackg round was gener- 



ated with MadGraph 4.4 461] for parton level processes, Pythia 6.4 47|] for hadronization 
and PGS-4 j^] for detector simulation. The total R parity-odd SUSY production cross 
section (citotai) for the low mass gaugino models are composed, to a first approximation, 
of only three contributions: production of chargino and the second lightest neutralino (i.e. 
cr^i^o/cTtotai; 47% ±2.5%); gluino pair production (i.e. cr^^/atotai; 28% ±3.3%); and chargino 
pair production (i.e. cr^±^T /crtotai; 23% ± 1.3%). The three sparticles produced with the 
largest production modes, namely g, xf^ ^ind X25 then decay with the dominant branching 
ratios shown in Table [III The ranges shown are for the subset of 700 models. The total 
SUSY production cross section is relatively large for this class of models given the relatively 
light gluino, charginos and neutralinos (o"totai = 9.65 pb±1.43 pb) over the set of 700 models. 
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Branching Ratios of the Low Mass Gaugino models in the Higgs-pole region 



Br{g X) % 


Br(x^ ^ X) % 


Br(xf ^ X) % 


UiUiX2 2 X (5.1 ±0.38) 
didiX2 2 X (5.0 ± 0.3) 
bbx2 15.1 ± 2.47 
UidiXi + h.c. 4 X (10.1 ± 0.75) 
tbxi +ii.c. 2 X (5.5 ±1.2) 


UiUiXi 2 X (12.5 ±0.57) 
didiXi 2 X (16.3 ±0.88) 
bbxl 16.1 ± 1.88 
/+rxi 3 X (2.9 ±0.49) 
i^iuiXi 3 X (5.7 ± 1.09) 


UidiXi 2 X (33.5 ±0.12) 
luixl 3 X (11.0 ± 0.07) 



TABLE II: Typical size of dominant branching ratios of the sparticles with the largest production 
modes emerging from proton-proton collision at the LHC over a subset of 700 models. Here li, d 
includes the first 2 generations of quarks and I includes all 3 generations of leptons (hence the 
factors of 2 and 3 in the Table). The factor of 4 includes u,d and the conjugate modes for the 
charginos. In addition to the three dominant sparticles arising from proton-proton collisions (the 
three cases considered in the Table), a small subset of models are found to produce light stops 
{rrii^ ~ 350 GeV) at the LHC which decay via ii — )• {ty^,bxi ^iXo) ^respectively, depending on the 
particular model point. 



The rather small variances around the central values for production cross sections 
and branching fractions suggest that the models in the Higgs-pole region are strikingly 
similar in their features, at least in terms of the phenomenology associated with the 
gaugino sector. This is not unexpected given previous studies of sparticle mass hierarchical 



patterns [l5, ll8|, |49|. As we will demonstrate in what follows, these similarities 



allow predictions to be made if excesses over SM background are observed at the LHC. 
Furthermore, as we will see in Section IIVI it is likely that these models will allow for a 
determination of the light gaugino masses and a partial determination of the neutralino 
LSP's eigencontent should a corroborating signal be observed in dark matter direct 
detection experiments. 



To illustrate the phenomenology of the low mass gaugino modelsin the Higgs-pole region 
we have chosen four benchmark models as presented in Table IHIi For each of these models 
we will compute the event rates for eight supersymmetric discovery channels defined by the 
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Key Spectra of Sample Models 



Til 

Label 


mo 


mi/2 


A 

Ao 


tan p 


rrig ruh 


Al 




rUg mi^ 


ruA — 


1 


2990 


148 


2503 


26 


476 119 


60 


117 


2959 1668 


2608 


2 


1238 


132 


-2007 


7 


407 116 


55 


109 


1250 421 


1467 


3 


2463 


133 


-2003 


50 


447 118 


58 


117 


2443 1353 


423 


4 


2839 


131 


-2401 


50 


451 119 


58 


118 


2812 1562 


355 



TABLE III: Four benchmarks to illustrate collider and dark matter signals of the low mass gaugino 
models in the Higgs-pole region. All models give a suitable relic density consistent with WMAP. 
Masses and dimensionful input parameters are given in units of GeV. The first and second gen- 
eration squarks are denoted by q. The top pole mass is set to 173.1 GeV and the sign of ^ is 
positive. Number in the table are rounded to the nearest integer. All values are computed with 
MicrOMEGAS 2.4 and SuSpect. 



following sets of cuts [24, |25| 



CUTCi : n{^) = 0, pt(ji) > 150 GeV, Pt{J2,J3,J4) > 40 GeV 
CUT Ca : n{e) = 0, n(6-jets) > 1 

n{j) > 4, ptUi) > 100 GeV, pr (j2, is, 74) > 40 GeV, > 0.2m,s 
n{j)> 4, pt{ji) > 100 GeV, m^s > 500 GeV 
n{j) + n{e) > 4, pt{3i) > 100 GeV, H^^ + $t> 500 GeV 
CUT Ce : = 3, n{]) > 2, pr {32) > 40 GeV 
CUTC7 : n{i) = 1, pT (Ji,j2,j3,j4) > 40 GeV, $t > 0.2meff 
CUT Cg : Z - veto, n{£+) = 1, n{e^) = 1, > 20 GeV . 



CUT C3 
CUT C4 
CUT C5 



(19) 



All eight channels involve a cut on transverse sphericity of St > 0.2 and a missing transverse 
energy cut of > 100 GeV, except for CUT Ci for which we impose > 150 GeV. 
Leptons of the first two generations (e, fi) are denoted collectively by £ and the number 
of leptons and the number of jets in an event are denoted by n{i) and n{j) respectively. 
Similarly, pri^i) and priji) refer to the transverse momentum of the i^^ hardest lepton 
or jet, respectively. The notation PriJi, j2, js, ji) means that the first through the fourth 
hardest jets in an event each have to individually pass the cut, and does not imply a sum. 
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LHC Significance for Channel d with 35 ph'^ and 1 fb-^ @ ^ = 7 TeV 





Jets iV|usY/\/^SM 


Leptons + Jets A^susy/V^sm 


Label 


CUT Ci 


CUT C2 


CUT C-i 


CUT d 


CUT Cs 


CUT Ce 


CUT C7 


CUT Cg 


1 


(2) [12] 


(1) [6] 


(2) [9] 


(2) [11] 


(2) [11] 


(0) [1] 


(1) [3] 


(0) [2] 


2 


(4) [21] 


(3) [14] 


(4) [21] 


(4) [24] 


(4) [23] 


(0) [2] 


(1) [6] 


(0) [1] 


3 


(3) [13] 


(1) [10] 


(2) [13] 


(3) [15] 


(3) [15] 


(0) [2] 


(1) [5] 


(0) [2] 


4 


(2) [15] 


(2) [10] 


(2) [13] 


(3) [16] 


(3) [15] 


(1) [2] 


(1) [5] 


(0) [2] 



TABLE IV: A^susy/\/^^ models of Table |TTI] for both (35 pb^i) and [1 fb^^] of integrated 

luminosity at the LHC with = 7 TeV. The (0) in the table means a significance of less than 1. 
We expect the entire set of our models discussed in Table |T] to surpass the 5(t significance threshold 
in jet-based channels early at LHC-7 with about an inverse femtobarn of data. 



If no value is specified for an object then no cut has been made for that object. In the 
specification of the cut Cs, the subscripts a and h indicate that the two opposite sign leptons 
may be of different flavors; a Z-veto is imposed on the invariant mass of the two leptons 
only in the case when they are of the same flavor, so as to avoid contamination from the Z 
boson peak produced through Standard Model production modes. 
We define the effective mass jries and H'^'' by 

4 4 

^eff = J^PT {ji) + $T, H^T = 5ZPT(Xi) , (20) 

i=l 1=1 

where Xj is a visible object (jet or lepton) and the summation, in both cases, is done over 
the first four hardest objects. The variable H^rp' is closely related to other definitions of Ht 



see 



50] for different definitions of Ht)- We define a model to be discoverable in a given 
channel (or for a given cut), Cj, if Ai'gusv ^ {"^^/N^, 10} , where -/Vgygy is the number 
of SUSY events and N^-^ is the number of background events. Further, we loosely refer to 
a 5(7 excess as one which satisfies iVgygy > Sa/A^, and a lower bound of ten events is 
imposed in rare cases where the SM background is insignificant for a specific channel. 

In Table IIVI we give an analysis of a broad range of event rates for the low mass gaugino 
models in the Higgs-pole region at y/s = 7 TeV with both 35 pb~^ and 1 fb~^ of luminosity 
under the cuts Ci as defined in Eq. f|T9|) . None of the models reach the discovery limit for 
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FIG. 2: (color online) Effective mass variable nics for the SUSY signal plus background with cut 
Ci at ^/s = 7 TeV. The SM background alone is shown shaded for comparison. For benchmark 1 
(top panel) , with a gluino mass of 476 GeV, we see a peak at rUeS = (725 it 25) GeV corresponding 
to a mass ratio of rn^^^^/mg = 1.52 it 0.055. For benchmark 2 (bottom panel), with a gluino mass 
of 407 GeV, a peak is observed at rries = (675 it 25) GeV which corresponds to a mass ratio of 
mlg^/nig = 1.66 ± 0.065. 
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the case of 35 pb~^. Benchmark point 2 has the largest significance for two reasons: It has 
the hghtest gluino mass of the benchmarks and the 2°*^ generation squarks are just above 
the TeV scale. Indeed, these models will produce discoverable signals with an increase of 
about a factor of 5 in luminosity, which may be expected within the next 6 to 8 months of 
data taking. However, any type of serious mass reconstruction will require about an inverse 
femtobarn of data. 

We find that the models analyzed produce a significant amount of jet events. These events 
arise from gluino decays via off shell squarks into fermion pairs with a chargino or neutralino, 
that is, g — > QiQixf and g — t- qiqiX2 with secondary 3-body decays X2 $t + 2 fermions 
and — )■ + 2 fermions. Additionally, one has a significant cross section for the 
direct production of charginos and neutralinos which can also give leptonic final states. 
Our analysis finds that the distribution of the transverse momentum of the hardest lepton is 
peaked near Pt(^i) = 20 GeV and falls off quickly near 60 GeV before imposing the cuts in 
Eq. ( fT9|) . The relatively soft leptonic decay products makes it more difficult to use leptonic 
signatures as discovery channels with limited data, as exhibited in Table IIVI However, the 
lepton + jets signal can be strong (see channel C5) where a large significance is achieved. 
Trileptonic signal Cg is only at the level of ~ 2a but would become visible with an increase 
in luminosity by a factor of six. The above features are generic to all models in the in the 
sample, given the rigid properties of the gaugino sector shown in Table HI 

The strongest signal of new physics will be in the multijet channel. In Figure Ej we plot 
the distribution in meg for two of our benchmark points using the cut Ci of Eq. ( |T9l) . The 
heavy solid line gives the supersymmetric signal events plus the SM background while the 
shaded area is the SM background. The peaks in this distribution can be identified with a 
typical accuracy of 25 GeV, which is half the bin size. A more statistically rigorous approach 
gives similar results. 



Several previous works [511] have shown that there is a relationship between the effective 
mass peak and the minimum mass of the gluino and the first two generation squark masses. 
Since in the low mass gaugino models that lie in the Higgs-pole region, the first and the 
second generation squark masses are always heavier than the gluino mass, the peak of the 
effective mass gives a relationship to the gluino mass. Analyzing the effective mass peak for 
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FIG. 3: (color online) Left: Distribution of the ratio of the effective mass peak to the gluino mass. 
The models plotted here are the 700 model subset and the peak is found after adding the SM 
background and applying cut Ci. We find the peak to be at 1.57 ± 0.085. Right: Distribution 
of the ratio of effective mass peak to the mass difference between the two lightest neutralinos 
under the same cut. The mass difference between the two lightest neutralinos corresponds to the 
upper bound of the edge in the OSSF dilepton invariant mass plot. We find the peak to be at 
12.50 ±0.721. 



cut Ci for all 700 simulated models we find in general 



P^^'^^l.Sm^ CUT Ci, (21) 



with the precise range being m^^^^/mg = 1.57 ± 0.085, as can be seen from the distribution 
in the left panel of Figure [3l We note that both of the benchmark cases in Figure [2] show 
this result explicitly. Thus a measurement of rn^s^^ provides an important early clue to the 
size of the gluino mass. Next, defining 

Am = m^o — m^o = (^a^o — 1 j m^o , (22) 

the mass relations found in Eq. ( fT3l) or Eq. ( fT4l) suggest that under cut Ci the peak in the 
effective mass distribution will be proportional to Am 

peak 

~ 1.5 X ^- = 1.5 X ^ — - . (23) 

Am {a^o - l)m^o {a^o - 1) 

The distribution of m^^^^/Am is shown to be peaked in the right panel of Figure [3l a result 
which follows from the left panel of Figure [3]and from the distribution in ag shown previously 
in Figure [TJ 
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OSSF dilepton Invariant Mass (GeV) OSSF dilepton Invariant Mass (GeV) 



FIG. 4: (color online) OSSF dilepton invariant mass for the SUSY signal plus SM background 
using cut C5 at ^/s = 7 TeV. The SM background is shown separately for comparisons. For the 
benchmark 1 (left panel) we see an edge at fn'^^- = 60 it 5 GeV and for the benchmark 2 (right 
panel) we see an edge at = 55 it 5 GeV, which agree well with the mass differences between 

the two lightest neutralinos in both cases, which are predicted to be 60 GeV and 55 GeV from 
theory (see Table [ml). 



The mass ratio plotted in the right panel in Figure [3] is noteworthy in that the quantity 
Am is measurable from the edge of the opposite-sign, same-flavor (OSSF) dilepton invariant 
mass distribution, ra^^- (for a recent study see 52]). In Figure H] we plot this distribution 



for the same two benchmark models from Figure [2] after applying the cuts C5 from Eq. (fT9l) . 
Upon reconstruction of the dilepton invariant mass for the two sample models, one observes 
clean edges near 55 GeV and 60 GeV for the two cases. For the complete set of the 700 
simulated models one finds 

/ \ I 0.75 m^o minimum 

edge ^ ^ ^ ( ^, ^ \ ^ — ) ^1 



< m^o - m^o = («^o - 1 J m^o = <j ' (24) 

1 . 7 m^o maximum 

In addition, from Eq. ( |T3l) we expect the upper bound of the OSSF dilepton plot to be less 
then 65 GeV which is the upper limit on Am found in the analysis which can be understood 
by using the appropriate predictions for the for each model point. 

In addition, because rrf^^- < Am, we can express the effective mass peak in terms of the 
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edge approximately as 

J 2 Ci^o — 1 , 

mS- ^ 3 X ^^Kt ■ (25) 

Thus we arrive at a very simple, but strong correlation between these two key observables 
at the LHC, i.e., m^+^'l and m^g^^- 

We therefore come to the conclusion that the low mass gaugino models in the Higgs-pole 
region are fully testable with early LHC data. If the models studied in this paper do indeed 
describe the supersymmetric content of our Universe, then the following three observations 
must follow: 

1. The dilepton invariant mass edge with an upper bound of (^(y^o — 1^ m^o < 65 GeV 
must be found. 

2. The multi-jet effective mass must be found, which peaks in the range 
550 GeV < mPg^"" < 800 GeV consistent with Eq.(^. 



3. The mass relation in Eq. fl25|) must hold. 
We now emphasize 

• LHC measurements can be used to estimate the dark matter mass in this model class. 
The upper bound of the edge in the OSSF dilepton invariant mass allows us to estimate 
the neutralino mass splitting and the scaling relation of Eq. allows us to infer the 
dark matter mass. 

• The effective mass peak and the dilepton invariant mass edge are strongly correlated 
via Eq. (2^) and provide cross-checks of the model. 

In the next section we will look for further avenues to exploit the remarkable predictivity 
of this model paradigm. 

IV. DARK MATTER DIRECT DETECTION EXPERIMENTS AND CONNEC- 
TION TO THE LHC 

The complementarity between dark matter detection experiments and collider signatures 



has been emphasized in many previous works (for a recent review see 



Here we will 



focus on this complementarity within the context of the low mass gaugino models in the 
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Higgs-pole region. We will show that experiments for the direct detection of dark matter 
such as XENON put further constraints on the parameter space of the model. 

We begin by noting that the predictions of Eqs. (I10fl2l) and the relic density con- 
straint largely ensure that the models yield predictions in narrow corridors as exhibited 
in Table HI Nevertheless, the properties of the neutralino, and in particular its scatter- 
ing cross section on nucleons, will depend on parameters such as fi, tan /3 and the re- 
sultant components riij which govern the wavefunction of the LSP. The features of the 
spin-independent neutrahno-nucleon scattering are easily understood in the models as they 
arise for large mo with the s-channel squark exchange suppressed and the scattering is 
dominated by Higgs exchange through the t-channel. Thus the spin independent scat- 
tering off target nucleus T arising via the interaction CiXiXiQiQi^ the limit of small 



momentum transfer is well approximated by o"|o'p 



fp/n XL 



given in 32 



53 



(p/n) 



Q '"p/ti _|_ _2_ jip/n) ^ ^ "m-p/^ 



27 J TG 



^q=c,b,t 



(4/i|or/7r)(Z/p +{A- Z)U)\ with 
C„'^ with the form factors /jf /^f^/") 



55( 1 and with coupling given by |53l-l55| 



92mg 
AmwB 



(^i[fl'2'^i2 - gynii]) DC 



1 1 



m 



H 



mi 



\ 2+ ' 

\mj^ m 



(26) 



The parameters 61^2 depend on eigen components of the LSP wave function and B, C, D 
depend on VEVs of the Higgs fields and the neutral Higgs mixing parameter a. For up 
quarks one has (5i, 62, B, C, D) = (riis, ni4, s^, Sa, Cq) and for down quarks {61, 62, B,C,D) = 
(ni4, — ni3, C/3, Ca, — Sq,). These simple relations reproduce numerical results of [32| and 
closely match the numerical work we do in this paper. 

For the four benchmark models of Table UTTl we present the spin- independent cross section 
of neutralino scattering on protons in Table |Vl However, from a survey over the collection 
of all the models in the Higgs-pole region we find a very broad range of possible scattering 
cross sections 

1' < < 4 X 10-^2 cm^ (27) 



4 X 10"^^ cm- 



XiP 



The largest of these are already ruled out experimentally from the null results of the 



CDMS II and XENON 100 experiments 30|, iMi- For the purposes of this paper we will 
assume a hard limit of a% < 6 x 10~^^cm^ for all neutralino masses under consideration 

XiP — 

as indicated by the XENON 100 experiment; this value is extremely conservative as their 
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Label 



Dark Matter and the Sample Models 



cjgpCm^ mi (B) ni2 {W) riis {Hi) nu (#2) ^CBuh?' 



1 


1.4 xlO" 


-46 


0.995 


-0.023 


0.093 


-0.015 


0.110 


2 


1.7 xlO- 


-46 


0.998 


-0.029 


0.058 


-0.012 


0.108 


3 


1.8 xlO" 


-44 


0.996 


-0.018 


0.092 


-0.012 


0.104 


4 


3.0 xlO" 


-44 


0.996 


-0.016 


0.085 


-0.011 


0.125 



TABLE V: Spin-independent cross section for neutralino scattering on protons for the benchmark 
models of Table IIIII Also given is the computed thermal relic density and the components nij of 
the LSP wavefunction. 

reported bounds are a factor of two more stringent, but we wish to allow for some uncertainty. 
A large fraction of the remaining models will be probed after longer exposures with XENON, 
or in future at other experiments. The distribution of our 12,000 models in the (m ^o, <^%j) 



plane is given in Figure |5] with both the CDMS II and XENON 100 limits indicated |30i |31|. 
Models which are being constrained by the XENON and CDMS data are those with 50 < 
tan/3 < 60. Note that the models in Figure [5] satisfy all the constraints discussed in Sec.fllTl). 

An important point to note is that dark matter direct detection experiments can be used 
to learn about soft supersymmetry breaking parameters. Figure [5] shows that once the spin 
independent cross section and neutralino mass are known from direct detection experiments, 
then mi/2 can be determined directly. Let us assume that a dark matter direct detection 
experiment observes a signal in the near future which is compatible with a neutralino LSP 
in the mass range 50 GeV ^ m^o ^ 65 GeV. Within the constraints of the of the Higgs-pole 
region even a crude measurement of the scattering cross section yields important information 
about the parameters of the model. The results shown in Figure [5] already demonstrate a 
correlation between cr|op and mi/2- For example a simultaneous estimation of m^o ~ 55 GeV 
and cr|op ~ 2 X 10""^^ cm^ would predict 125 GeV ^ mi/2 ^ 140 GeV due to the correlated 
nature of the parameters within the Higgs-pole region. This, in turn, would have testable 
consequences for the gaugino sector at the LHC. 

The XENON bound can be mapped into a constraint on tua- This constraint is more 
restrictive than the one from collider bounds. Without direct detection constraints, a pseu- 
doscalar mass as low as itla — 190 GeV is allowed, as it satisfies the Tevatron search limits as 
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FIG. 5: (color online) The spin independent cross section cr|op versus neutralino mass. Points are 
colored according to the value of mi/2 taken. Applying the XENON and CDMS limits we see that 
mi/2 is preferred in the 120 GeV to 155 GeV region. 



well as the indirect constraints imposed above. For example, one such model in Figure |5] has 
rriA = 190 GeV, tan /3 = 56, m^o = 60 GeV, riu = 0.994 and riia = 0.102; for this particular 
model, cr?o ~ 5.5 x 10"^^ cm^ in excess of what is allowed by XENON 100 data. Thus the 
XENON constraint is stronger than the Tevatron bound for this point. More generally, we 
obtain a limit arising from the dark matter direct detection constraint: 

ruA Z 300 GeV XENON Constraint . (28) 

Including uncertainties in the form factors that enter the computation of cr?Op one may 
loosen or tighten this constraint a bit; however, the point here is that the constraints on ttla 
become rather strong from the XENON data. The value quoted above is particular to the 
requirements within the confines of the scaling predictions in Eq. f lTU]) and the mass range 
Eq. ( JT2I) . However, other models with radiative electroweak symmetry breaking are also 
strongly constrained. We have performed a separate analysis to investigate minimal super- 
gravity models which satisfy the WMAP constraints of Eq. (ITTj) via stau-co-annihilation, 
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XENON Dark Matter Search, Tevatron and LHC 




CP Odd Higgs Mass GeV 



FIG. 6: Displayed is the small subset of the 12,000 models which are those corresponding to large 
tan /3 and with low niA within reach of LHC-7 in the first year (a majority of the 12,000 models 
have heavier niA and lie off this graph). The LHC estimated projected reach (magenta curve) with 



ISO. 



ated tau pairs and 6-tagging 



63] is indicated. Models ruled out by the XENON 100 experiment 



30] are in blue (squares) and we have taken a conservative cut cj^Op < 6 x 10^ pb to account 
for theoretical and experimental uncertainties. Red (diamonds) are allowed models in this mass 
range of (myi,tan/3). Constraints on sparticle mass limits as well as other constraints are imposed 
as discussed in Section m however the models ruled out by the Bg — )• n^fj.^ constraint are shown 
explicitly in green (circles) to illustrate its effects. The shaded yellow region indicates where the 
Tevatron has excluded niA- We conclude that the XENON 100 constraints are very severe in this 
part of the parameter space. 



which have a heavier neutralino mass than the models studied here (owing to mass limits on 
the stau) and we find that the present XENON data imposes only a slightly weaker lower 
)ound of niA ~ 250 GeV. Constraints of this type have been studied in SUGRA models 



561 ] and in generic weak scale MSSM models in references |57l . l58l ] and more recently in the 



context of low mass dark matter in references 59l-l61| . The results presented here show that 
for dark matter in the 50 GeV region, the constraints on the CP-odd Higgs sector in models 
of radiative breaking are also quite strong. We anticipate that the lower bound on will 
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only get stronger as additional data from XENON arrives (for projections see e.g. 62|). 

It is interesting to note that Eq. fl28l) is precisely the mass scale for which the LHC will 
be sensitive to the production of the pseudoscalar Higgs with 1 fb^^ at = 7 TeV 63[. It 
is therefore possible to probe the pseudoscalar Higgs at LHC-? in the 2r + 6 — tagged jets 
channel within a subset of the models. In conjunction with the measurements of Eq. (l25l) this 
could serve to extract the value of tan /3. We therefore exhibit the subset of the 12,000 models 
with large tan /3 in Figure |6] and plot tan /3 vs. the CP-odd Higgs m^. The heavy black line 
(yellow shaded region) is the Tevatron direct search limit, while green points are eliminated 
from Tevatron constraints on Bt{Bs — fi^fi^). Blue squares represent models that are 
eliminated by the (conservative) imposition of cr|op < 6 x 10"^'* cm^ from XENON 100 
results. The red points are the surviving models with < 550 GeV and the estimated 
LHC-7 reach for 1 fb~^ is shown by the (solid) nearly vertical magenta curve. We note that 
there are a number of cases which could give detectable signals at the LHC, and in addition, 
a substantial portion these models correspond to spectrum with a light CP-odd Higgs mass 
which have a neutralino mass and spin independent cross section that lie close to the range 
of observation relevant to the XENON experiment. 



V. DETERMINING GAUGING AND HIGGSINO CONTENT OF THE LSP AND 
THE SOFT PARAMETERS FROM THE INTERSECTION OF DARK MATTER 
AND LHC DATA 

In this section we will further connect the LHC to dark matter detection. In particular 
the data from both types of experiments can be combined to extract information on the soft 
SUSY breaking parameters as well as the eigencontent of the neutralino LSP. 



A. Determining Eigencontent of the Neutralino LSP 

Let us assume that dark matter direct detection experiments have determined (or at 
least restricted) the possible range of LSP mass and spin-independent cross section. Unfor- 
tunately, in the models, this information leaves the LSP eigencontent in terms of gaugino 
and Higgsino components still undetermined. The model points in Figure |5]in the {m^o,a% ) 
plane that are unconstrained by the XENON data have large fluctuations in their Higgsino 
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FIG. 7: (color online) Higgsino eigencontent of the LSP displayed as a function of a^o = m^o/m^o. 
The models are indicated by the gluino mass. Once a^o is measured via corroborating evidence at 
the LHC and in dark matter detection, and the gluino mass is is deduced at the LHC, the Higgsino 
eigencontent, y^|ni3p + |ni4p may be determined. 

content. Hence we need to turn to LHC data in concert with direct detection data in order 
to sort this out. The two types of measurements at the LHC required are the ones discussed 
above: a measurement of the edge in the OSSF dilepton invariant mass and a measure- 
ment of ml^^^. Taken together with dark matter detection results, these quantities can help 
determine the eigencontent of the LSP as we now show. 

Previously, we have seen that a measurement of the edge in the OSSF dilepton invariant 
mass at the LHC gives us an upper bound on Am, the mass difference between the two 
lightest neutralinos (see Eq. (!24|) ). Taken together with the LSP mass measured by dark 
matter experiments as well as the LHC, this information then gives an experimental deter- 
mination of the mass ratio m^o/m^o = 1 + Am/m^o. This quantity is the horizontal axis in 
Figure [71 

Additionally, a measurement of fn^^^^ at the LHC gives us a good estimate of rrig, as 
can be seen in Eqs. f l?T|TB]) . In Figure [7] we have shaded the model points according to 
the value of rrig. Hence, given this information together with the value of m^o/m^o along 
the horizontal axis allows us to estimate the Higgsino fraction of the LSP plotted along 
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FIG. 8: (color online) cj?0p versus m-o/m^o distributed in Ao/mo- The ratio Aq/itiq exists in 
separate regions relative to m^o/m^o and knowledge of the OSSF edge at the LHC can point to 
the soft parameter space. 

the vertical axis. Thus one can then essentially read off the Higgsino eigencontent of the 
neutralino dark matter from Figure [3 Clearly this determination will be rough due to 
uncertainties at every stage, but it provides a first step in the determination of the gaugino 
vs. Higgsino eigencontent of the LSP. 

In complementary fashion, once dark matter experiments can measure the lightest neu- 
tralino mass one can then determine Og as well (see Eq. [25]) . Finally, we note that in the 
limiting case when the models approach the pure bino limit for the neutralino, it is seen 
from Figure [7] that the ratio of the second lightest neutralino to the LSP approaches 2 and 
the gluino mass is driven towards its lowest value. In summary, these observables combined 
together would lend strong support for the model class. 

B. Determination of Ao/mo 

The density of possible values of the ratio of soft SUSY breaking parameters Ao/itlq in 
the models from our scan is shown in Figure [HI on a plot of versus m^o/m^o . Let us 
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now assume that dark matter experiments have determined a^Op as well as m^o. One can 
see that current bounds on a^o^ as discussed in the previous section already rule out some 
ranges of Aq/itio. Most of the remaining models congregate around Aq/ttio ~ ±1. 

As in the previous subsection, the horizontal axis m^o/m^o can be determined by a 
combination of LSP mass obtained from dark matter experiments (as well as LHC) together 
with Am determined from a measurement of the edge in the OSSF dilepton invariant mass 
at the LHC. As future bounds on improve, some further information on Ao/ttiq will be 
attained. 

It is interesting that the combination of the two types of experiments could help determine 
the scalar trilinear Aq relative to rriQ as the trilincar couplings arc otherwise difficult to 
measure from the LHC data alone. As an explicit example to the above general statements, 
models with Ao/mo — 1 are found to have a mass splitting between X2 and xl of 45 GeV < 
Am < 50 GeV . The models with Am near the upper hmit of the allowed range favors 
the opposite case with Ao/mo ~ — 1. In addition, majority of the models which congregate 
around a^o — m^o/m^o < 1.8 are ruled out by XENON. Indeed, as emphasized in the 
previous section, the LHC should be able to determine the dark matter mass of any of the 
models with the largest uncertainty at about the 20% level. 



VI. CONCLUSION 



We have analyzed a predictive model relevant to early SUSY discovery at the LHC at 
y/s — 7 TeV. We claim that within the framework of minimal supergravity unification, 
models with ~ 50 GeV dark matter must be found in early LHC data, or they will be ruled 
out. Our analysis was targeted at the mass scale where the LSP can have a mass of this size 
consistent with astrophysical and particle physics constraints, and where the relic density of 
dark matter is largely governed by the presence of the light CP even Higgs pole. Connected 
are the mass of the relic lightest neutralino, and the gluino mass, the latter of which has an 
upper bound of about 575 GeV in this model class. Such a gluino can be detected in the 
early runs at the LHC from its distinctive decay signatures consisting of energetic leptons 
and jets along with a sizeable missing energy. 

The model can be further checked in direct detection experiments such as XENON via a 
detection of event rates consistent with the spin independent neutralino- proton cross section 
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(T?Op which has a theoretical upper bound near 10 cm^ while a large collection of these 
models tend to be in the range a% = 10~^^^^ cm^. In connection with the above, we 

Xi'P 

showed that the current experimental limits from XENON 100 already put limits on the 
model and lead to a lower bound on the CP-odd Higgs mass of > 300 GeV, which is 
more stringent than the current constraints from direct searches for the production of the 
pseudoscalar from the Tevatron. 

It was further shown that measurements of certain signatures at the LHC can allow one 
to estimate the neutralino mass and the gluino mass with the LHC data. With sufficient 
luminosity the kinematic edge in the OSSF dilepton invariant mass distribution directly 
allows one to estimate the neutralino dark matter mass due to scaling in the gaugino sector; 
namely the ratio of the masses of two lightest neutralinos are related by a scale factor, and 
this scale factor is close to 2. Similarly, from the distribution, one can infer the gluino 
mass. 

If the low mass gaugino models within the Higgs-pole region studied in this paper do 
indeed describe the supersymmetric content of our Universe, then there are three absolute 
predictions which must be found in the data. First, the location of the dilepton invariant 
mass should be seen in a narrow range near 50 GeV. Since this mass edge is very close to the 
mass of the dark matter particle, its measurement will determine the dark matter mass to 
~ 20%. Second, the multijet effective mass under our cuts will peak in the range 550 GeV ^ 
m^^^^ ^ 800 GeV. Since this peak is related to the gluino mass via m^^'^^ ~ 1.5 x nig, this 
measurement will give a first estimate in the determination of the mass of the gluino. Third, 
we have deduced a simple relation between the peak in the effective mass and the dilepton 
invariant mass edge via Eq. f l25|) that can be checked directly with LHC data. 

In addition, it was shown that the intersection of constraints from the LHC and direct 
detection experiments provide further information about the SUSY model. A combination 
of accelerator and direct detection data sets can provide estimates of tan /3 and Ao/rrio; 
can tell us about the gaugino and Higgsino content of the dark matter; and can provide 
information about the mass of the dark matter particle. The model class is consistent with 
the very recent ATLAS and CMS data with 35 pb^^. A most exciting feature of the analy- 
sis given here is that the required data to test the model will be taken in the very near future. 
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